Topologically Protected Zero Modes in Twisted Bilayer Graphene 
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We show that the twisted graphene bilayer can reveal unusual topological properties at low ener- 
gies, as a consequence of a Dirac-point splitting. These features rely on a symmetry analysis of the 
electron hopping between the two layers of graphene and we derive a simplified effective low-energy 
Hamiltonian which captures the essential topological properties of twisted bilayer graphene. The 
corresponding Landau levels peculiarly reveal a degenerate zero-energy mode which cannot be lifted 
by strong magnetic fields. 



PACS numbers: 73.43.Nq, 71.10.Pm, 73.20.Qt 

I. INTRODUCTION 

One of the most fascinating aspects of graphene is its 
band structure which can be fundamentally changed in 
several different ways by modifying its lattice structure. 
This happens because the honeycomb lattice of mono- 
layer graphene has two independent sublattices and the 
electron, as it moves through the lattice, has to change its 
sublattice and hence the character of its wavefunction^ 
Thus, even small local changes in the lattice structure 
lead to the appearance of gauge potentials which are as- 
sociated with the phase of the electronic wavefunction 
in each sublattice. As a consequence, there is a one-to- 
one correspondence between graphene's structure with 
the topological features of the electronic states. 

Another amazing property of graphene is its honey- 
comb structure which yields an electronic low-energy ef- 
fective theory which is Lorentz-invariant in two dimen- 
sions (2D), and thus corresponds to 2D Dirac fermions. 
This Lorentz invariance is robust because the energy 
associated with sublattice coupling, that is, the inter- 
sublattice hopping energy t (w 3 eV) is the dominant 
energy scale in the system. Lorentz invariance persists 
even when the lattice is modified either by external forces 
(strain, shear, etc.) pi by external fields, or by the ad- 
dition of more layersP For instance, AB-stacked bilayer 
graphene is described, at low energies by two sets of mas- 
sive Lorentz invariant Dirac particles per valley and spin. 
In the simplest models where only nearest neighbor hop- 
pings are taking into account, there is still an accidental 
degeneracy that makes the particle-like band of one fla- 
vor to be degenerate with the anti-particle- like (hole-like) 
band of the other flavor at the K (K') point in each valley. 
This degeneracy can be easily lifted by the application of 
a perpendicular electric field that breaks the inversion 
symmetry in the system^ Although the Lorentz invari- 
ance is preserved, the wavefunction of the electrons at 
low energies is modified— whereas in monolayer graphene 
the Dirac fermions carry a Berry phase ±7r, the Berry 
phase is ±27r in AB-stacked bilayer grapheneP At low en- 



ergies, trigonal warping splits this "double" Dirac point 
into three with a Berry phase 7r and an additional one 
with — 7rp"l a situation that persists for a translational 
mismatch between the layersP 

Twisted bilayer graphene, in which the two layers have 
a rotational mismatch described by an angle as com- 
pared to the perfect AB stacking, is another example 
where lattice structure and wavefunction topology are 
directly interconnected. In fact, from the experimen- 
tal point of view, twisted graphene is more the rule 
than the ex ception. It naturally occurs at the surface of 
graphitep^l in graphene grown in the surface of SiCP or 
graphene grown by chemical vapor deposition on metal 
substrates^ As compared to monolayer graphene, each 
Dirac flavor is then split into two copies that are sep- 
arated in reciprocal space by a wave vector AK, as a 
function of 0.^1 Inter-layer ho pping results in a renor- 
malization of the Fermi velocit y 10 " 13 ! as well as in a van 
Hove singularity at relativel y low energy as compared to 
that in monolayer graphene! 9 ! 10 ! Once again, Lorentz in- 
variance is preserved at low energies but the nature of the 
electronic wavefunctions is modified in a profound way. 

In this paper, we investigate the topological aspects 
of the band structure of twisted graphene bilayers at low 
energies in the continuum approximation, for small angle 
mismatches as compared to perfect AB stacking. We 
identify possible topological classes that describe band 
inversion symmetry. These topological classes determine 
the relative Berry phase between the two copies of Dirac 
particles, which have either the same or opposite Berry 
phase. If the two Dirac cones are related by time-reversal 
symmetry, such as the K and K' points in monolayer 
graphene, the Berry phases are naturally opposite. In 
twisted graphene bilayers, however, the two Dirac points 
emanating from different layers are not time-related, and 
the symmetry of the inter-layer hopping term enforces the 
Berry phases to be identical. We show that this feature 
yields a topologically protected zero-energy Landau level, 
in contrast to the former case. The scenario may be 
tested in quantum Hall measurements. 
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FIG. 1: First Brillouin zone for twisted bilayer graphene. The 
1BZ of the upper layer (dashed hexagon) is rotated by an 
angle 6 with respect to that of the lower layer (full hexagon). 
The corners, where Dirac points occur, are labeled by Kg 
and K^'\ respectively. 



The paper is organized as follows. In Sec. [TTJ we dis- 
cuss the model of twisted bilayer graphene and the sym- 
metry properties of their bands (Sec. II A) that are fixed 



by the form of the inter-layer hopping. Furthermore, we 
present an effective two-band model (Sec. II B ) that dis- 



plays the same topological low-energy properties as the 



original four-band model. Section [TTT] is devoted to the 
discussion of the Landau-level spectrum in twisted bi- 
layer graphene, in the perspective of possible quantum- 
Hall measurements. 



II. MODEL OF TWISTED BILAYER 
GRAPHENE 

If we neglect, for the moment, hopping between atoms 
in different layers, the electronic properties of twisted bi- 
layer graphene are described by two copies of the Hamil- 
tonian for monolayer graphene (we use units with % = 1): 



H (k) = V F 



k* 
k 



(1) 



iky is the 2D 



where vf is the Fermi velocity and k = k, 
wave-vector relative to the K (K') points of the rotated 
layers (see Fig[T]). For a twist (rotation) angle 9 ^ 0, 
each of the two inequivalent Dirac points, which reside 
at the corners of the first BZ K and K' , are split into two, 
separated by a wave vector AK = K Kg, where Kvm 
is the position of the K point in the lower (upper) layer 
and — K(e) the position of the points K' and K' e , respec- 
tively. Throughout this paper, we will work in the con- 
tinuum limit around a single pair (K,Kg) and, therefore, 
neglect commensuration effects between the two layers 
which could enlarge the unit cell in position space and 
thus fold it back in reciprocal space. This procedure is 
mostly justified because the coupling between the pairs 
of Dirac conef^EH is negligible.^ 

The Hamiltonian describing the electronic properties 
of twisted bilayer graphene reads 

/i/ (k+AK/2) H ± \ 

H W-[ H [ ff (k-AK/2)J' (2) 



where H± is the hopping matrix between the two layers. 
Equation ^ refers to an expansion around the Q point of 
Fig. [T] The analysis of the Moire pattern formed by the 
twisted bilayer shows that, for small angle 9, the hopping 
matrix H± may have three different fo rms c orresponding 
to the three main Fourier components! 10 * 15 ^ This leads to 
three different types of inter-layer hopping terms, 



#X = *-L I 1 1 J > H ± = *-L [ e ±i<t> e T*< 



(3) 



where <fi — 2tt/3 and t ± is a hopping parameter which 
generally depends on ^(MUEU 



A. Symmetry of the bands and Berry phases 

In contrast to a lattice Hamiltonian that may be ana- 
lyzed with the help of global symmetries, such as time re- 
versal or lattice inversion, H(k) in Eq. (|]) is a continuum 
model in which the two Dirac points are no longer related 
by these symmetries. However, H(k) and the inter-layer 
hopping term H± may be investigated via symmetries 
that involve directly the energy bands, such as rotation, 
mirror, and inversion symmetry. Whereas for H± = 
the rotation and inversion symmetries are respected, the 
latter are broken for non-zero inter-layer hopping, and 
we therefore restrict the discussion to the inversion sym- 
metry I of the bands. We emphasize that this inversion 
symmetry is defined with respect to the bands in recip- 
rocal space, in contrast to a previous analysispUn which 
the more common definition of inversion symmetry with 
respect to the lattice was used. 

In the absence of any inter-layer hopping, the Hamilto- 
nian |2]) would be ambiguous since one could work equally 
well with Hq in the second layer. As a consequence, there 
remain two possible representations of inversion (with 
different spinorial expressions) Ii and 1% : 

2~i : H(-k) = -H(k) or 2~ 2 : H(-k) = -H(k)*. 

The minus signs in both transformations maps positive to 
negative energy states at opposite wave vectors and vice- 
versa: -E'(k) — > — E(— k), whereas the complex conjuga- 
tion in I2 changes the relative phase between the spinor 
components and, hence, the Berry phase of a cone. 

The phases of the two Dirac points are thus different 
for the two representations. For the I\ case, the Berry 
phases at fixed energy are opposite such that a merging 
transition of the two Dirac points (at 9 = 0) corresponds 
to a vanishing total Berry phase and consequently to the 
possible opening of a band gap. This situation arises 
e.g. in the framework of the model discussed in Refs. 
1181 which describes monolayer graphene under strong 
strainP In contrast to this rather well-known topologi- 
cal universality class, the transformation I2 yields Berry 
phases that are the same at fixed energy E and opposite 
to those at —E and thus represents a second universality 
class for Hamiltonians describing pairs of Dirac points. 
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The transformation I2 may be represented in terms of a 
tensor product of Pauli matrices, X 2 = 0® ® cr^, where 
vti,, and Cw,, describe the intra-layer and the inter-layer 
spinorial spaces, respectively. The phases are eventually 
fixed by the symmetry of the interlayer hopping H± , and 
the particular forms ([3| happen to be invariant under X^. 

Notice, however, that as a consequence of the inter- 
layer hopping terms ([3]), the two Dirac points are no 
longer exactly at the same energy, but we neglect this en- 
ergy shift here because it is associated with a very small 
energy scale ~ 1 meV. As a second-order perturbation, 
this indeed scales like ^/(■upAK) with t± of the order 
of 100 meV and v F AK ~ 1 eVP 



B. Effective two-band model 



Energy (a.u.) 




FIG. 2: (Color online) Generic band structure for the low- 
est four bands of the twisted bilayer graphene, around the 
Q point, with only Hj_ considered. The inset of the fi gur e 
pictures the effective bands of the model Hamiltonian (13 1. 



Both band structures are inversion-invariant with respect to 
the Q point. 



In the opposite limit, t± ^> v F AK, the model (2l may 
be reduced to an effective two-band model, similarly to 
the case of a perfectly AB-stacked (6 — 0) graphene 
bilayer P This is done in two steps. First, one replaces 
H±_ in Eq. ([2| by a simplified inter-layer hopping term, 



H 



err 



t± 





1 



(5) 



in the limit where t± 3> v F AK, i.e. for small tilt an- 
gles. The inter-layer hopping term ^ is reminiscent of 
the Bernal bilayer case. In spite of this simplification and 
the difference in the energy scales, the resulting Hamil- 
tonian may be viewed as a representative of the topolog- 
ical universality class that also includes the original four- 
band model. We consider the eigenvectors of Eq. ^ 
in terms of the 4-spinor basis {ipA,'4 ) B,'4 ) A l ,ipB'} where 
(A, B) belongs to the two sublattices of the first layer and 
(A 1 , B') to that of the second one. With the particular 
form of Eq. ([5]), the zero-energy sector is spanned ipA 
and ipB't whereas tpB and ipA' are strongly hybridized 
by the inter-layer hopping term. Their symmetric and 
anti-symmetric combinations are the eigenstates at — t± 
and t±, respectively. In order to describe the electronic 
properties in the vicinity of E = 0, one may therefore 
project the Hamiltonian onto the reduced {ipA,ipB'} ba- 
sis and neglecting terms of the form Eijj B / A ,, which are a 
product of the energy E ~ and the small components 
ipB I A 1 ■ The eigenvalue problem reads 

v F (k + AK/2)*ip B = E^ A (6) 

v F (k + AK/2)tp A +i±^A' = Evjj B ~0 (7) 

hipB + v F (k - AK/2)*i> B , = Et/j a > (8) 

v F {k-AK/2)ip A , = Eijj B i . (9) 

By rewriting Eqs. Q and 

-v F /t ± (k + AK/2)ij A = ^ A > (10) 
-v F /t ± (k-AK/2)*i> B , = i> B , (11) 



and substituting them into Eqs. (|6j) and ([£]), one obtains 
the Schrodinger equation 



if eff (k)' ^ 



4>B> 



= E 



IpA 
IpB' 



(12) 



in terms of the effective two-band Hamiltonian 



H ctt (k) = -J 



t ± V k 2 - (AK/2Y 



(AK*/2) 2 




(13) 

which is similar to that for the nematic transition of the 
interacting Bernal graphene bilayerP^l 

The generic form of the band structure obtained from 
the effective two-band model ( 13 ) is depicted in the in- 
set of Fig. [2] One notices the two Dirac points origi- 
nally situated at K and Kg, separated by a wave vector 
I AK| ~ 6 /am. the first BZ, in terms of the intra-layer dis- 
tance a — 0.142 nm between neighboring carbon atoms. 
In order to see the linearity of the dispersion relation in 
the vicinity of these contact points, one can further ex- 
pand the Hamiltonian (5) around ±AK/2, by defining 
k = q ± AK/2, with |q| < |AK|/2. This expansion 
yields two Dirac Hamiltonians 



± 



IAK 



q* 
q 



(14) 



with identical chirality for the two contact points. 

Furthermore, the bands have saddle points at k = 
between the two Dirac points. The effective two-band 
Hamiltonian therefore captures the logarithmic van-Hove 
singularity in agreement with previous theoretical^ and 
experimental studies.^ 

The most salient features of the Hamiltonian ( 13 ) 
are its chiral properties. In agreement with the above- 
mentioned general symmetry considerations, electrons at 
a fixed energy at the K-point have the same chirality 
as those at the Kg point, such that in both cases the 
electron experiences the same Berry phase 7 = it (and 
7 = —it at the points K' and K e ) on a closed orbit 
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around one of the Dirac points. This is obvious in the 



E(eV) 



low-energy expansion (14) around the two Dirac points. 



As for AB-stacked bilayer graphene, the Berry phase ac- 
quired on an orbit enclosing both Dirac points is then 
27, as one may also see from the limit of vanishing twist 
angle (9 = 0) which reproduces the perfectly AB-stacked 
bilayerPFurthermore, it becomes apparent from the form 



of the Hamiltonian ( 13 1 that the merging of the Dirac 



points is not accompanied with a gap opening, in contrast 
to the Hamiltonian discussed in Ref. 18, which desribes 
the same band structure in the semi-metallic phase but 
which belongs to another topological class, described by 
the symmetry I\ . 



ILL LANDAU LEVELS OF TWISTED BILAYER 
GRAPHENE 

One of the most prominent consequences of these topo- 
logical properties is the presence of a (doubly-degenerate) 
zero mode that emerges in the presence of a quantizing 
magnetic field, in which case the Hamiltonian ( 13 1 may 
be written in terms of the usual ladder operators a and 
a) , with 



.t] = 1: 



H, 



a^ 2 — a 



a — a 
2 



(15) 



where up = 2vpeB /t± is the cyclotron frequency and 
a = AKIb/^VZ- As compared to the perfectly AB- 
stacked bilayer (a = 0), where one readily obtains the 
Landau level (LL) spectrum^ one notices that the addi- 
tional terms couple states only of the same parity. One 
obtains thus two classes of eigenstates tp2n and ip2n+i 
that may be written in the usual harmonic-oscillator ba- 



sis 



|m), with a^a\m) = m\m): 



2n(+l) 



E ( i m{+1) ) |2m(+l)>, (16) 

where the components § x m and § 2 m are to be determined 
recursively. 



A. Zero-energy levels 



Before discussing the LL spectrum of Hamiltonian 



(15), we investigate the zero-energy states, which may 



be obtained analytically from the equation Hb^ = 0. 
As for the AB-stacked bilayer, one obtains two distinct 
solutions, one with even and one with odd parity, such 
that the zero-energy level is orbitally two-fold degenerate, 
in addition to its usual four-fold spin-valley degeneracy 
and the orbital degeneracy described by the flux density 
riB = eB/h. The two zero-energy states, which are remi- 
niscent of coherent states, are directly obtained from the 
secular equation: 



*o = A/" cosh(a*aT) 







0) 



(17) 
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FIG. 3: (Color online) Landau levels of twisted bil aye r 
graphene, obtained from the numerical solution of Eq. (15 1. 



The characteristic energy scale of the van-Hove singularity has 
been chosen as |q| 2 ojc = (AA") 2 /8m = O.leV. The red line 
indicates the zero-energy level, and the numbers correspond 
to the filling factors in the gaps. 



and 



sinh(a*aT) 





|n = 0) 



(18) 



in terms of the normalization factors Ao/i- These states 
are generalizations of the zero-energy states (0, |0)) and 
(0, |1)) of the AB-stacked bilayer P. We also notice that 
we neglect any other potential lifts due to Zeeman effect 
and/or interactions P2 The existence of those zero-mode 
states is independent of the strength of the magnetic field. 
This is to be compared to the other topological class of 
Dirac cones with opposite Berry phases where the de- 
generacy of the zero-mode is lifted.^ This protection is 
solely determined by the topology of the model Hamil- 
tonian since the band-structures are otherwise identical. 
Notice further that the scaling of the Landau levels rela- 
tive to AK is different for the two topological classes. 



B. Full Landau-level spectrum 

The full LL spectrum, which has been calculated nu- 
merically, is depicted in Fig. [3j In the regime of small 
magnetic fields, the LLs with small index n <C nc, 
where nc denotes roughly the LL which crosses the 
van-Hove singularity, \Jn c B ~ vp\ AK\ 2 /4:^/2tj_, display 
the V Bn behavior, which is the benchmark of massless 
Dirac fermions, as expected for the linear dispersion be- 
low the van-Hove singularity. Because of the two Dirac 
points, these LLs, as well as the zero-energy level dis- 
cussed above, are twofold degenerate, in addition to the 
fourfold spin-valley degeneracy. The twofold degeneracy 
due to the Dirac-point splitting in the twisted bilayer is 
lifted once the LLs approach and eventually cross the 
van-Hove singularity, n > nc, above which the LLs scale 
as ~ B(n + 1/2), as one would expect from the parabolic 
dispersion relation revealed by Hamiltonian ( 13 1 at high 
energies. 



Notice, however, that as for AB-stacked bi- 
layer graphene the effective two-band approximation is 
no longer valid at higher energies (far beyond the van- 
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Hove singularity) because of the presence of the remain- 
ing bands, which become visible there. 

The LL spectrum in Fig. [3] allows one to understand 
the main features of a quantum Hall effect (QHE) in 
twisted bilayer graphene. The topologically protected 
zero-energy LL, with its altogether eightfold degeneracy, 
yields a QHE at filling factors v = ±4 (taking into ac- 
count spin degeneracy). This feature is independent of 
the interlayer hopping strength t± and of the van-Hove 
singularity, which is triggered by the twist angle 6. For 
other LLs, the position of the van-Hove singularity deter- 
mines their degeneracy If the LLs remain well below the 
singularity, n -C nc, they maintain their eightfold de- 
generacy, and one would therefore expect Hall plateaus 
at filling factors v = ±4(2n + 1) = ±4, ±12, ±20, but 
one would expect additional plateaus at v = ±8, ±16, ... 
for LLs with n > nc- The experiment!^ indicates, even 
at rather low magnetic helds of B ~ 5T, only the zero- 
energy LL is eightfold degenerate, whereas a plateau at 
v = ±8 has been observed. This stipulates that the en- 
ergy of the van-Hove singularity is as small as the first 
excited LL. 



IV. CONCLUSIONS 

In conclusion, we have investigated the topological 
band structure of twisted bilayer graphene, in the frame- 



work of a symmetry analysis of the inter-layer hopping 
in the continuum limit. For small and moderate twist 
angles 8, the two copies of the Dirac point (that are not 
related by time-reversal symmetry) are described by the 
same Berry phase, due to the the symmetry of the inter- 
layer hopping term. Therefore, they belong to a different 
topological class than the (usual) two Dirac points, which 
are related by time-reversal symmetry. In the presence 
of a quantizing magnetic field, these particular topolog- 
ical properties yield a protected zero-energy LL with an 
eight-fold degeneracy that may be evidenced in quantum 
Hall transport measurements. 
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